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Abstract 

Geometrical structure of matter at subnuclear densities is investigated in the 
presence of a degenerate gas of neutrinos as encountered in stellar collapse. The 
crystalline phases with spherical, cylindrical and planar nuclei as well as with spher- 
ical and cylindrical nuclear bubbles are considered by using a compressible liquid- 
drop model. This model allows for uncertainties in the lepton fraction 1l in addition 
to those in the nuclear surface tension -Bsurf and in the proton chemical potential 
in bulk neutron matter fip\ The phase diagrams obtained at zero temperature 
show that only the phases with rod-like and slab-like nuclei occur at typical values 
of Yl; Esuri and fip \ whereas the bubble phases, especially with spherical bub- 
bles, are at best expected at hypothetically low values of Yl and/or -Bgurf- For the 
rod-like and slab-like nuclei, thermally induced displacements are calculated from 
their respective elastic constants. It is found that at temperatures appropriate to 
supernova cores, thermal fluctuations would destroy the layered lattice of slab-like 
nuclei almost independently of the nuclear models and of the degree of the neutrino 
degeneracy. 

PACS : 26.50.+X; 97.60.Bw 

Keywords : Dense matter; Ground state; Thermal fluctuations; Stellar collapse 
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1 Introduction 



To clarify the mechanism of collapse-driven supernovae (see, e.g., Refs. [0, ^ for review) 
necessitates understanding of the properties of matter under extreme conditions char- 
acterized by high densities and temperatures as well as of its interactions with electron 
neutrinos, which are produced by electron capture processes on the nuclei present during 
the collapse of an iron core of a massive star. Diffusion of the neutrinos thus produced is 
controlled by their coherent scattering off the nuclei via neutral-current weak interactions 
0, 1^. It is generally accepted that at densities above ~ 10^^ g cm~'^, where the dynamical 
time scale of the collapse is small compared with the diffusion time scale, these neutrinos 
are trapped, forming an ultrarelativistic and degenerate Fermi gas. This plays a role in 
establishing (3 equilibrium in the systemQ and, until the core ceases to collapse and re- 
bounds, in increasing the system temperature T up to order or greater than 10 MeV. The 
density region in which such a trapping occurs and the lepton fraction, i.e., the number 
of the leptons of electron flavor present per nucleon locally in this region: 



where rib, and riy are the nucleon, electron and electron- neutrino number densities, 

respectively, depend on the rates of the electron capture and of the neutrino scattering. 

These rates are in turn affected by changes in the states of supernova matter such as 

changes in the nuclear composition and a melting transition from the phase with the 

nuclei to the phase of uniform nuclear matter. 

^ Hereafter, the corresponding (3 equilibrated material is referred to as supernova matter. 
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At subnuclear densities and at temperatures below the critical temperature of order 10 
MeV, nuclear matter is separated into a high- density liquid part, usually characterized by 
roughly spherical nuclei, and a low-density gas part, composed mostly of neutrons. When 
the density is close to the saturation density of symmetric nuclear matter, ps ^ 2.7 x 10^^ 
g cm~^, the system is expected to have rod-like and slab- like nuclei as well as rod-like and 
roughly spherical bubbles, as originally indicated by Ravenhall et al. and Hashimoto 
et al. using a liquid drop model in the zero-temperature approximation. These earliest 
investigations show that such non-spherical nuclei and bubbles, forming a Coulomb lattice, 
arise from a delicate balance between the interfacial energy and the electrostatic energy 
including the lattice energy. With increasing density, a bcc lattice of spherical nuclei 
changes into a two-dimensional triangular lattice of rod-like nuclei and then into a one- 
dimensional layered lattice of slab-like nuclei, until the liquid part and the gas part begin to 
be replaced by each other. After further transformations to a two-dimensional triangular 
lattice of cylindrical bubbles and then to a bcc lattice of spherical bubbles, the system 
turns into uniform nuclear matter. This behavior was confirmed by Lassaut et al. |^ in 
their three-dimensional Thomas-Fermi calculations at finite Yl and T relevant to stellar 
collapse. Those kinds of non-spherical nuclei and of bubbles are often referred to as 
nuclear "pasta" since they are similar in shape to spaghetti, lasagna, etc. 

The presence of non-spherical nuclei and of bubbles is expected to have consequence 
to the hydrodynamics and neutrino transport in supernova cores; these are ingredients of 
a simulation study of stellar collapse, bounce and explosion (for recent progress in this 
study, see, e.g., Ref. 0). Let us first note that the neutrino wavelengths, typically of 
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order 20 fm, are comparable to or even greater than the internuclear spacing, leading to 
diffractive effects on the neutrino elastic scattering off such a periodic spatial structure 
of nuclear matter |^. These effects, induced by the internuclear Coulombic correlations, 
would act to reduce the scattering rates and hence the lepton fraction Y^. For the bcc 
lattice of spherical nuclei, such a reduction was examined by Horowitz [P] by calculating 



m 



the associated static structure factor. As clarified by Bruenn and Mezzacappa [TO 
their numerical calculations of stellar collapse, however, the resultant window that allows 
low-energy neutrinos to diffuse away from the core during its collapse is too narrow to 
drastically suppress the degree of the neutrino degeneracy. It is also noteworthy that 
non-spherical nuclei and bubbles are elongated in specific direction. In such direction, the 
neutrino scattering processes are no longer coherent, in contrast to the case of roughly 
spherical nuclei whose finiteness in any direction yields constructive interference in the 
scattering. The last point to be mentioned is that the changes in the nuclear shape are 
accompanied by discontinuities in the adiabatic index, denoting how hard the equation of 
state of the material is. As Lassaut et al. |0] suggested, these discontinuities may influence 
the core hydrodynamics during the initial phase of the collapse. 

For neutron star matter, which has vanishing neutrino degeneracy, we examined the 
density region in which the phases with non-spherical nuclei and with bubbles are preferred 
over the bcc phase of spherical nuclei and the homogeneous phase by generalizing a 
compressible liquid-drop model advanced by Baym, Bethe and Pethick |jT^ (hereafter 
BBP) so as to incorporate uncertainties in the nuclear surface tension -Esurf and the 
proton chemical potential /i^^*^ in the dripped neutron gas |1T[. We found that as -Egurf 
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decreases or fx^^ increases, such a density region is enlarged appreciably. Fortunately, 
those uncertainties are of less importance to supernova matter at subnuclear densities, 
which will be studied in the present paper. This is because after the onset of the neutrino 
trapping, the resultant neutrino degeneracy keeps the denser part of nuclear matter stable 
against further neutronization; not only does the proton fraction x in this part remain 
typically of order 0.3, not far from the regime accessible to experiment, but also the 
dripped nucleon gas, arising mainly from thermal effects, is too low for the values of fx^^ 
to make any essential difference in the density region occupied by nuclear "pasta". To 
be taken into account in analyzing supernova matter, in addition to the nuclear model 
dependence of the stability of the "pasta" phases, is its dependence on the trapped lepton 
fraction Fl- The values of 11, depending on the sizes and shapes of the denser part which 
produces, scatters and absorbs electron neutrinos via weak interactions, exert feedback 
on the condition for the stability of the "pasta" phases. Note, for example, that with 
increasing 11, the melting density pm approaches ps- This is due to the tendency ||T^ 
that as the denser part is at larger neutron excess, the saturation density is lowered. By 
allowing for uncertainties in Yl coming partly from such a feedback effect, we set a range 
of Fl as 0.2 < 1l < 0.4, rather wider than the usually adopted values ranging 0.3-0.38 (see 
Refs. fl], 0] and references therein). The equilibrium phase diagrams of zero-temperature 
supernova matter are drawn for this range of 1l; we thus find that at typical values of 
Esuri, Pp"-* and Yl, the phases with rod-like and slab-like nuclei lie between the usual bcc 
phase and the uniform phase, as in the case of neutron star matter ||1 1|| . 

For these two "pasta" phases, we evaluate thermally induced displacements of the 
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nuclei from their elastic constants, which were described by Pethick and Potekhin |]H 
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terms of the interfacial and electrostatic energies. They noted that the low-dimensional 
spatial order exhibited by the phases with cylindrical and planar nuclei is analogous to 
that of the liquid crystals, i.e., columnar phases and smectics A, respectively. According 



to our previous analysis of neutron star matter fill], thermal fluctuations may melt the 
layered lattice of planar nuclei rather than the triangular lattice of cylindrical nuclei at 
temperatures appropriate to matter in neutron star crusts. In the case of supernova 
cores of interest here, as we shall see, not only the temperature scale but also the elastic 
constants are considerably large as compared with the case of neutron star crusts. The 
temperatures and elastic constants thus increased play a role in enhancing and suppressing 
the fluctuational displacements, respectively. The resultant rivalry between these two roles 
is finally examined. 

In Section 2, a compressible liquid-drop model for nuclei and bubbles, together with 
the associated equilibrium conditions, is designed to describe zero-temperature matter at 
subnuclear densities and various lepton fractions. In Section 3 we therefrom derive the 
equilibrium phase diagrams for changes in the nuclear conformation. Properties of such 
conformation changes are discussed in Section 4. In Section 5 we estimate displacements 
of rod-hke and slab-like nuclei at finite temperatures. Section 6 is devoted to conclusions. 
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2 Model for supernova matter at subnuclear densi- 
ties 



In this section, we construct the free energy of supernova matter for densities near the 
melting density pm, and write down the conditions for its equihbrium. We concentrate 
on the influence of the lepton fraction 1l on the free energy. Since the finite-temperature 
effects on the free energy, i.e., nuclear and electronic excitations, nucleon evaporation 
from the denser part of nuclear matter, multiple shapes and sizes of this part, etc., are of 
less importance, the temperature T is taken to be zero. All we bear in mind is thus to 



combine the BBP-type compressible liquid-drop model used by us [|11[] for the description 
of zero-temperature, inhomogeneous nuclear matter embedded in a sea of electrons with 
the energy contribution of electron neutrinos. At T = 0, the less dense part of nuclear 
matter is composed entirely of the neutrons dripped out of the denser part. 



Following Ref. we consider five phases characterized by the shapes of the denser 
part of nuclear matter, i.e., sphere, cylinder, slab, cylindrical hole and spherical hole, 
respectively; each phase is composed of a single species of nucleus or bubble at a given 
baryon density n\, and lepton fraction Yl- The Wigner-Seitz approximation is adopted in 
describing the energy of the respective Coulomb lattice: A cell in the bcc lattice, including 
a spherical nucleus or bubble of radius tn, is replaced by a spherical Wigner-Seitz cell 
having the same center and radius a cell in the two-dimensional triangular lattice, 
including a cylindrical nucleus or bubble having an infinitely long axis and a circular 
section of radius tn, is replaced by a cylindrical Wigner-Seitz cell having the same axis 
and a circular section of radius r,,; a cell in the one-dimensional layered lattice, including 
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a planar nucleus with width 2rN, is equivalent to a planar Wigner-Seitz cell having the 
same central plane and width Ir^- The values of Tc for these phases are chosen so that 
each Wigner-Seitz cell may have zero net charge. 

2.1 Energy of matter 

Let us now regard each species of leptons (nucleons) as uniformly distributed everywhere 
in the system (inside and outside the nuclei or bubbles), and write down the total energy 
density -Etot of the system as 



Here, wn is the energy of the nuclear matter region (the region containing protons) in a 
cell as divided by the cell volume; wl is the lattice energy per unit volume as given for 
the nuclei or bubbles having nonzero tn; is the number density of neutrons outside 
the nuclei or inside the bubbles; i?n, -Ee and Ey are the energy densities of the neutron 
matter, of the electron gas and of the neutrino gas, respectively; u is the volume fraction 
occupied by the nuclei or bubbles. 

In Eq. (2), the expressions for wn, wl, -^n and are given by Eqs. (3), (8), (10) and 
(11) in Ref. ||ri|. A new component, the neutrino energy density, is expressed in the form 
of an ideal ultrarelativistic Fermi gas: 



where ky = {Q'K'^ny)^ is the neutrino Fermi wave number. Recall that the energy of the 
nuclear matter region is written in the form of a compressible liquid-drop model as (Eq. 




(nuclei) , 



(bubbles) . 



(2) 




(3) 
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(3) in Ref. ^) 

un[{l — x)mn + xmp]c^ + unW{k, x) 

+Wsuriin, X, Tin, tn, u, d) + Wc{n, X, tn, u, d) (nuclei) , 



WN(n,a;,nn,rN,rc,rf) 



[1 — u)n[{\ — x)m^ + xmp]c^ + (1 ~ u)nW{k, x) 
+Wsuriin',x,nn,rt^,u,d)+wc{n,x,rf^,u,d) (bubbles) , 

(4) 



where mn {fnp) is the neutron (proton) rest mass; n is the nucleon number density in the 
nuclear matter region; W{k,x) is the energy per nucleon for uniform nuclear matter of 
nucleon Fermi wave number k = (37r^n/2)^/^ and proton fraction x; Wgurf is the nuclear 
surface energy per unit volume; wq is the self Coulomb energy (per unit volume) of protons 
contained in a cell; d is the dimensionality defined as c? = 1 for slabs, d = 2 ioi cylinders 
and (i = 3 for spheres. This expression for includes three parameters Ci, C2 and C3, 
associated with uncertainties in the proton chemical potential /j,^^ in pure neutron matter 
as contained in W{k,x) and those in the nuclear surface tension £'surf(= ^t^Wsurf/ud). Ci 
determines the magnitude of fi^^ (not including the rest mass) as (Eq. (4) in Ref. [|n|]) 

;.(°) = -C,nT . (5) 
We basically set Ci = 400 MeV fm^; this case is consistent with the generally accepted 



behavior among recent literature as exhibited in Fig. 1 of Ref. |Tl|. As the case may be, 
Ci = 300, 500, 600 MeV fm^ will be also taken for the sake of comparison. The other two 
parameters C2 and C3 are defined as (Eq. (6) in [pT| ) 



E,„rf = C2tanh(^-^jEf,?P , (6) 

where fx^J^^ = dE^/dn^ — m^c^ is the neutron chemical potential in the neutron gas (not 
including the rest mass), and EfJ^f is the BBP-type surface tension (Eq. (7) in Ref. WM)- 



As shown in Ref. [jTT| for the ground-state matter with 71^ = and Ci = 300, 400, 500, 600 
MeV fm^, to set C2 = 1-0 and C3 = 3.5 MeV well reproduces the Hartree-Fock (HF) 
results obtained by Ravenhall, Bennett and Pethick |15| (hereafter denoted by RBP) using 
a Skyrme interaction. It is important to examine whether or not this tendency persists 
for the ground-state supernova matter; in this case, the proton fraction x of interest is 
considerably large relative to the case of = in which it is confined to x < 0.2. When 
we still set C2 = 1.0 and C3 = 3.5 MeV, as can be seen from Fig. 1, the surface tension 
calculated from the equilibrium conditions to be described in Subsection 2.2 agrees fairly 
well with the RBP results for Ci = 400 MeV fm^ and 1l = 0.2, 0.3, 0.4.0 Such agreement 
has been confirmed also for Ci = 300, 500, 600 MeV fm^ and 1l = 0.2, 0.3, 0.4. In the 
present work, we thus give C3 a fixed value, 3.5 MeV, and C2 a range of values including 
unity, 0.01-3. 

2.2 Equilibrium conditions 

Zero-temperature supernova matter with nuclei or bubbles of given shape, in its equilib- 
rium, fulfills the conditions for stability of the nuclear matter region against change in the 
size, neutron drip, /5-decay and pressurization, as in the case of neutron star matter (see 
Section 2.2 in Ref. [|ll|])- These conditions come from minimization of the energy density 
Etot with respect to five variables n, x, nn, tn and u at fixed lepton fraction 1l given by 
Eq. (1) and baryon density rib given by 

_ J un + {1 — u)nn (nuclei) , 

~ 1 (1 - u)n + urin (bubbles) , ^ ' 



^At fixed Yl and for densities of interest here, the proton fraction x is Hmited to a rather small range 
exhibited in Fig. 1. 
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as well as under charge neutrality, 

^ _ f xnu (nuclei) , 

^ I xn{l — u) (bubbles) . 

The resultant expressions for the size, drip and pressure equilibria are the same as those 



given for neutron star matter (see Eqs. (14), (15) and (21) in Ref. [Tl[|). Only the (3- 
equilibrium condition (Eq. (18) in Ref. [jll|) is modified by the neutrino degeneracy; it 
now reads 

(J-e- ^u- (fTln - mp)c^ = /X^^^ - /xjf ^ , (9) 

where 

/ijy = hkiyC (10) 

is the neutrino chemical potential; /ie, fx^^ and fi^^ are the electron chemical potential, 
the neutron chemical potential in the nuclear matter region and the proton chemical 
potential in the nuclear matter region, respectively (see Eqs. (19), (16) and (20) in Ref. 



At given nb, 1l and nuclear shape, we have calculated the minimum value of Etot 
for electrically neutral supernova matter. Such calculations begin with elimination of /ij, 
from condition (9) with the help of Eqs. (1) and (10), leading to a cubic equation for fii . 
The solution for /Xe can then be written in terms of n, x, Uj^ and by using the Cardan's 
formula. For the rest of the process that allows us to find the equilibrium values of n, x, 
rin, and Vc and thus the minimum value of -Etot, we can follow a line of argument of 
Ref. Q. 

At the same nb and 11, the energy density of uniform nuclear matter in /5-equilibrium 
with leptons has also been evaluated. The total energy density of this uniform system is 

11 



written as 

Etot = nW{k, x) + n[(l - x)mn + xmp]c^ + E^{ne) + E^{n^) , (11) 

where the presence of muons is ignored. We then minimize Etot with respect to x at fixed 
baryon density nb{— n) and lepton fraction 11 as well as subject to charge neutrality, 

Tie = xn . (12) 

The resultant /3-equilibrium condition, given by 

dW(k,x) . X 2 

— ^ -iJ,e + iJ,u + {mn-mp)c , (13) 

leads to the optimal value of x and hence of i?tot- By comparing this value of £'tot with 
the values obtained for the five crystalline phases, we have determined the phase giving 
the smallest energy density at various values of ^b, 1l, Ci and C2, and thereby drawn the 
phase diagrams for the ground-state supernova matter, as plotted in Figs. 2 and 3. 

3 Phase diagrams 

Let us first examine the fi^^ and -Egurf dependence of the density region of each phase; 
such dependence can be typically seen from Fig. 2 that exhibits the phase diagram drawn 
at Fl = 0.3 and Ci = 300, 400, 600 MeV fm^ over the rib versus C2 plane. We have thus 
confirmed the feature as mentioned in Section 1 that the density region where neither 
the usual bcc phase nor the uniform phase is stable is almost independent of fj,^^ because 
of negligibly small nn. We remark in passing that the melting density pm increases with 
increasing Ci (lowering n^^), in contrast to the case of neutron star matter (see Fig. 
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3 in Ref. [|rT|]). This is because such a change in helps the nuclear matter region 
keep the protons within itself, rather than the neutron gas region share the protons. The 
-^'surf dependence, on the other hand, is essentially the same as that obtained for neutron 
star matter: Not only is the density region occupied by nuclear "pasta" reduced with 
increasing C2, but also the sequence of the nuclear shape for the lower -Esurf {C2 ~ 0.3), 
i.e., sphere — ^ cylinder slab cylindrical hole spherical hole uniform matter 
(with increasing n?,), changes into sphere cylinder — slab uniform matter for the 
typical -Esurf (C2 — 1-0). However, it is to be recalled that uncertainties accompanying C2 
are smaller for supernova matter than for neutron star matter (see Section 1). 

We proceed to consider the Yl dependence of the phase diagrams by observing Fig. 
3 drawn for Fl = 0.2,0.25,0.3,0.4 and for Ci = 400 MeV fm^. For these values of 
Yl, within the typical -Esurf region (C2 — 1.0), the phase with slab-like nuclei dissolves 
into uniform matter without undergoing a transition to the phase with cylindrical or 
spherical bubbles. The most viable situation in which at least one of the bubble phases 
is energetically favored for a finite range of density is expected in a rather hypothetical 
parameter region, Yl — 0.2 and C2 — 0.7. It is interesting to note that as Yl increases, 
the melting density pm goes up and becomes close to ps- As mentioned in Section 1, this 
is due to the proton fraction dependence of the saturation density; for details, see Fig. 1 
in Ref. [|l^ that shows the BBP model for the energy of bulk nuclear matter underlying 
W{k,x) used here. We also find that the density region occupied by nuclear "pasta" 
is considerably larger than that displayed by Fig. 3 in Ref. ||ll| for neutron star matter 
having a lepton fraction of order 0.03 and a surface tension of order 0.1 MeV fm^^. This 
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is because such a density region is controlled by the tendency towards reduction in the 
total surface area, which is in turn determined by -Esurf and hence by Yl (see Fig. 1). This 
tendency communicates with the fact that for supernova matter, the energy difference 
between two successive phases is generally of order 1-10 keV per nucleon (see Fig. 4), an 
order of magnitude larger than that for neutron star matter (see Fig. 4 in Ref. UTTH ). 

We conclude this section by asking why the bubble phases do not appear in the equilib- 
rium phase diagrams for our typical nuclear model (Ci = 400 MeV fm^ and C2 = 1.0) in 
contrast with the results obtained by Lassaut et al. [0 in their three-dimensional Thomas- 
Fermi calculations using a Skyrme nucleon-nucleon interaction. As expected from a simple 
liquid-drop argument of the electrostatic and interfacial effects ignoring the bulk energy 
[0, they found that all the phases with non-spherical nuclei and with bubbles occur for 
a finite range of pressure (or density). We ascribe this contrast primarily to the differ- 
ence in the adopted properties of bulk nuclear matter at a proton fraction of about 0.3; 
description of these properties relies on interpolation between the fairly-well-determined 
properties of pure neutron matter and symmetric nuclear matter. This ascription is sup- 
ported by the speculation that the role of curvature corrections in stabilizing the phases 
with bubbles relative to the phases with nuclei, which was allowed for in Ref. but has 
been disregarded in the present work, is not sufficient to cause the bubble phases to ap- 
pear in the vicinity of C2 = 1.0 on the phase diagrams. Such insufficiency can be viewed 
from Fig. 5, in which the curvature energy gain per nucleon for spherical bubbles over for 
planar nuclei and the difference in energy per nucleon between the phases with spherical 
bubbles and with planar nuclei have been compared for 1l = 0.3, Ci = 400 MeV fm^ and 
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6*2 = 1-0.0 The curvature energy for given conformation is expressed per nucleon as 

- ^ r ' (bubbles) , 

where Uc is the curvature thermodynamic potential per unit length. We have set the 
values of Uc, which have yet to be determined well even at proton fractions near x = 0.5, 
as ujc = 0.1,0.2,0.5 MeV fm~^ by reference to the estimates, ~ 0.1-0.2 MeV fm~^ at 



X ~ 0.3, made by Kolehmainen et al. in the Thomas-Fermi theory with Skyrme 



interactions. Figure 5 indicates that for the typical parameter set leading to x ~ 0.3, the 
value of Uc required by the presence of the bubble phases is at least ~ 0.5 MeV fm~\ 
considerably larger than the estimates referred to. For Uc ~ 0.1-0.2 MeV fm~^, all the 
curvature corrections contribute to the vicinity of C2 = 1.0 on the phase diagrams is to 
lower the transition density from sphere to cylinder and that from cylinder to slab by 
several tens percent and to accordingly enlarge the density region surrounded by the bcc 
phase and the uniform phase. 

4 Properties of phase transitions 

At a pressure at which the transition between one nuclear shape and another occurs, 
various quantities such as rih, Tc tn, n, x and are more or less discontinuous because 
of the dependence of the interfacial and electrostatic energies on the dimensionality d 



^ Note that for this parameter set and at densities close to the melting point pm, the phase with 
cylindrical bubbles is less favorable than that with spherical bubbles even in the absence of the curvature 
corrections which lower the energy for the spherical-bubble phase further than that for the cylindrical- 
bubble phase. 
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(see Eqs. (5) and (8) in Ref. [Tl|]). As can be derived from Fig. 4 by implementing the 
double-tangent constructions denoting the coexistence of the two neighbouring phases, 
the discontinuity in amounts to of order 1%. Due to such a weak first-order nature of 
the transition, the transition point, properly characterized by a constant pressure, is well 
described by a constant density determined by the comparison between the total energy 
densities for the respective phases. 

Let us then examine the changes in the sizes of the nucleus or bubble and of the 
Wigner-Seitz cell, tn and Vc, associated with the structural transitions. In Fig. 6 we have 
plotted the values of and evaluated as functions of n^, for 1l = 0.3, Ci = 400 MeV 
fm^ and C2 = 0.01, 1.0, 2.5. Discontinuities in tn and Tc at the transition points have been 
thus clarified. In the case of C2 = 0.01 in which all the five crystalline phases appear, 
we observe the same dimensionality dependence of as observed for neutron star matter 
(see Fig. 5 in Ref. [|11]), i.e., it is largest for li = 3 (spheres and cylindrical holes) and 
smallest for d = 1 (slabs). 

We finally address the question of how the neutron densities outside and inside the 
nuclear matter region change as the system melts into uniform nuclear matter. Compar- 
ison of these densities with the one for the uniform phase has been made in Fig. 7 for 
yL=0.3, Ci=400 MeV fm^ and C2 = 0.01,0.1,1.0,2.5. We can observe from this figure 
that for densities up to the melting point, the density of the dripped neutron gas remains 
negligibly small compared with the neutron density inside the nuclear matter region. This 
behavior, which stems from the neutrino trapping (see Section 1) and manifests itself for 
the other values of 1l and Ci, is consistent with the results obtained by Ogasawara and 
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Sato |T8| in their Thomas- Fermi calculations at various values of 1l and zero temperature. 
We also find from the comparison of Fig. 7 with Fig. 6 in Ref. that the difference in 



the neutron density profiles between the crystalline and uniform phases at the melting 
point is remarkable in contrast with the case of neutron star matter, in which the neu- 
tron density profile is fairly levelled just below the melting point. In supernova matter, in 
which neighbouring nuclei (or bubbles) appear to touch and fuse each other at the melting 
point, neutrino elastic scattering off nuclear "pasta" is thus expected to take effect until 
it becomes uniform. 

5 Thermal fluctuations 

Given occurrence of the phases with rod-like and slab-like nuclei as expected for typical 
values of -Esurf and Yj^ from the zero-temperature phase diagrams exhibited in the previous 
section, we turn to the estimates of thermally induced displacements of these nuclei from 
their equilibrium positions. In making such estimates, we utilize the expressions for such 



displacements that were written in Ref. |11| on the basis of the elastic constants |14 



controlling small deformation of the nuclei. For the layered lattice of slab-like nuclei, the 
mean square of the displacement v of layers in their normal direction reads (Eq. (36) in 
Ref. 0) 

(H')-T^S^lnf^^) , (15) 



where is the Boltzmann constant, a = 2rc is the layer spacing, B and Ki are the elastic 
constants given by Eqs. (28) and (29) in Ref. [11], respectively, and the length scale L of 



the lattice is assumed to be far larger than r^. For the two-dimensional triangular lattice 
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of rod-like nuclei, the mean square of the displacement vector v of cylinders in a plane 
perpendicular to their elongated direction is expressed as (Eq. (37) in Ref. |Tl|) 



I |2\ ^ 



where A = y2Ks/{B + 2C), B, C and K3 are the elastic constants given by Eqs. (31), 
(34) and (35) in Ref. |jlll, respectively, a = (27r/v^)^/^rc is the lattice constant for 
hexagonal cells, and the linear dimension of the lattice in the plane including v is assumed 
to be much larger than Tc and to be much smaller than the length of the nuclei. It 
is important to notice that all the elastic constants, derived by Pethick and Potekhin 



Tj] from the calculations of the electrostatic and interfacial energy increase due to the 
nuclear deformation, are proportional to the electrostatic energy density, wq+l, given by 
Eq. (8) in Ref. |Tl|. Since wq+l in turn behaves as oc x"^, the non-spherical nuclei in 



supernova cores are more difficult to deform than those in neutron star crusts. As was 
noted in |]TT| concerning the crossover temperature Tq from the quantum-fluctuation to the 



thermal-fluctuation regime, expressions (15) and (16) hold for T ^Tq but underestimate 
the displacements for T ^ Tq. In the case of supernova matter considered here, Tq can be 
roughly estimated to be of order 10^ K for planar nuclei and of order 10^ K for cylindrical 
nuclei. 



Let us then ask how large the root-mean-square displacement yd'^P) of a planar or 
cylindrical nucleus is relative to the shortest distance, a/2 — tn, between the surface of 
the nucleus in its equilibrium position and the boundary of the cell containing it. The 
ratio of these two lengths for the planar (cylindrical) nucleus has been plotted in Fig. 8 
(9) at T = 10 MeV and L = 10 km, appropriate to the subnuclear density regime inside 



supernova cores; the L dependence of the displacement, which is logarithmic, makes little 
differences. We can find from these figures the dominance of the relative displacement 
of the planar nucleus over that of the cylindrical nucleus, yielded by the logarithmic 
factor appearing in Eq. (15), as well as the suppression of the relative displacements 
by -Esurf, which were obtained also for neutron star crusts (see Figs. 7 and 8 in Ref. 
|1TT|). This -Esurf dependence then communicates with the suppression of the relative 
displacements by Yl as exhibited in Figs. 8 and 9 via the relation between -Esurf and Yl 
shown in Fig. 1. In examining the competition between the influences of temperature and 
of elasticity as pointed out in Section 1, it is noteworthy that the relative displacements 
thus obtained are larger than those obtained for neutron star crusts by less than a factor 
of ten. This is because the influence of the difference in the typical temperature by two 
orders-of-magnitude surpasses that of the above-mentioned difference in the elasticity of 



the nucleus. In case we set the melting condition as yd^P) = 0-/2 — tn, the enhanced 
relative displacements in supernova cores do not modify our previous conclusion |TT 



obtained for neutron star crusts that the triangular lattice of rod-like nuclei is stable 
against thermal fluctuations for our typical nuclear model, but change the conclusion 
about the layered lattice of slab-like nuclei from its marginal instability |[Tl| depending 
on the nuclear models to its instability almost independent of them. Given uncertainties 
in the temperature profile of the core, it is useful to estimate the critical temperature 
at which the relative displacements become unity; the results for Tc have been plotted in 
Figs. 10 and 11 for the layered lattice and the triangular lattice, respectively. 
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6 Conclusion 



The dependence on the lepton fraction Fl of the phase diagrams associated with the 
changes in the nuclear shape has been studied for zero-temperature supernova matter at 
subnuclear densities. We have found that with increasing 11, the density region occupied 
by the "pasta" phases, i.e., the phases with non-spherical nuclei and with bubbles, becomes 
larger. Despite the considerable increase in the lepton fraction from neutron star matter, 
leading to enhancement of the proton fraction in the nuclear matter region and suppression 
of the density of the dripped neutron gas, not only the sequence of the nuclear shape 
(sphere — > cylinder — > slab — > uniform with increasing density) obtained for the typical 
nuclear model but also the effectiveness of thermal fluctuations at the dissolution of slab- 
like nuclei has been shown to be analogous to the case of neutron star matter. 

The actual presence of nuclear "pasta" in supernova cores depends on how the un- 
certain properties of bulk nuclear matter having a proton fraction of order 0.3 affect the 
stability of the "pasta" and on whether or not the time needed to form the "pasta" is 
small compared with the time scale of stellar collapse. If future work about these ques- 
tions were in favor of the occurrence of the "pasta" , it would be significant to calculate the 
finite-temperature equation of state and neutrino transport in the presence of the "pasta" 
phases and to incorporate them into a supernova simulation with a neutrino radiation hy- 
drodynamics code. This would clarify the influence of the "pasta" on the lepton fraction 
at the onset of the neutrino trapping. 
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Fig. 1. The surface energy per unit area (the surface tension) as a function of x. 
The thick hues are the present results obtained from Eq. (6) for Ci = 400 MeV fm^, 
C2 = 0.01,0.1,0.5,1.0,2.5 and = 0.2,0.3,0.4. The sohd curve is the RBP resuh from 
their HF calculation |jl5|, and the dotted curve is the BBP result HTsf. The cross denotes 
the value adopted by Ravenhall et al. p. 

Fig. 2. Zero-temperature phase diagrams on the rib versus C2 plane, evaluated for 

= 0.3 and Ci = 300, 400, 600 MeV fm^. 

Fig. 3. Zero-temperature phase diagrams on the ny, versus C2 plane, evaluated for 
= 0.2, 0.25, 0.3, 0.4 and Ci = 400 MeV fm^. 

Fig. 4. Energy per nucleon relative to that of uniform matter calculated for Yl = 0.3, 
Ci = 400 MeV fm^ and C2 = 1.0 as a function of baryon density n^. The symbols SP, 
C, S, CH and SPH stand for sphere, cylinder, slab, cylindrical hole and spherical hole, 
respectively. 

Fig. 5. Energy per nucleon of the phase with spherical bubbles relative to that with 
slab-hke nuclei, calculated for Yl = 0.3, Ci = 400 MeV fm^ and C2 = 1.0 as a function of 
baryon density nb. The absolute value of the curvature energy per nucleon of the phase 
with spherical bubbles is also evaluated for the curvature thermodynamic potential per 
unit length uj^ = 0.1, 0.2, 0.5 MeV fm'^ 

Fig. 6. Size of a nucleus or bubble, tn, and of a Wigner-Seitz cell, Tc, calculated for 
Yl = 0.3, Ci = 400 MeV fm^ and C2 = 0.01, 1.0, 2.5 as a function of baryon density rib- 
The symbols SP, C, S, CH and SPH stand for sphere, cyhnder, slab, cylindrical hole and 
spherical hole, respectively. 
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Fig. 7. The neutron densities obtained for Fl = 0.3, Ci = 400 MeV fm^ and C2 — 
0.01,0.1,1.0,2.5asa function of baryon density ^b- The hues classified by [n(l — x)]n and 
n-a represent the neutron densities in the nuclear matter region and in the neutron gas 
region for the phases with nuclei, respectively. The lines classified by [n(l — a;)]uni denote 
the neutron densities in uniform nuclear matter. 

Fig. 8. The root-mean-square displacement of a planar nucleus at k-eT — 10 MeV, 
divided by the shortest distance between the surface of the nucleus in its equilibrium 
position and the boundary of the cell containing it. The curves are obtained for Ci = 400 
MeV fm^ and for various sets of 1l and C2 as a function of baryon density n^. The thick 
curves lying between the two vertical lines are the results in the density region in which 
the phase with planar nuclei is energetically stable. 

Fig. 9. The root-mean-square displacement of a cylindrical nucleus at k^T — 10 MeV, 
divided by the distance between the surface of the nucleus in its equilibrium position and 
the boundary of the cell containing it. The curves are obtained for Ci = 400 MeV fm^ 
and for various sets of 1l and C2 as a function of baryon density n^,. The thick curves 
lying between the two vertical lines are the results in the density region in which the phase 
with cylindrical nuclei is energetically stable. 

Fig. 10. The critical temperature for the phase with planar nuclei as a function of 
baryon density nb. The thick curves lying between the two vertical lines are the results 
in the density region in which the phase with planar nuclei is energetically stable. Ci is 
taken to be 400 MeV fm^. 
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Fig. 11. The critical temperature for the phase with cyhndrical nuclei as a function 
of baryon density rib- The thick curves lying between the two vertical lines are the results 
in the density region in which the phase with cylindrical nuclei is energetically stable. Ci 
is taken to be 400 MeV fm^. 
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